Effects of longer-range interactions on unconventional superconductivity 
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We analyze the effect of the non- vanishing range of electron-electron repulsion on the mechanism of 
unconventional superconductivity. We present asymptotically exact weak-coupling results for dilute 
electrons in the continuum and for the 2D extended Hubbard model, as well as density-matrix 
renormalization group results for the two- leg extended Hubbard model at intermediate couplings, 
and approximate results for the case of realistically screened Coulomb interactions. We show that 
Tc is generally suppressed in some pairing channels as longer range interactions increase in strength, 
but superconductivity is not destroyed. Our results confirm that electron-electron interaction can 
lead to unconventional superconductivity under physically realistic circumstances. 
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I. INTRODUCTION 

That unconventional superconductivity can arise in 
models with short-ranged repulsive interactions between 
fermions has been known for some time, beginning with 
the pioneering work by Kohn and Luttingei— . This issue 
has been revisited multiple times in the last few decades, 
since the discovery of unconventional superconductivity 
in the cupratcs. For the Hubbard model with local re- 
pulsion U and fcrmionic bandwidth W, the existence of 
superconductivity in p, d, f or g— wave, as well as sign- 
changing s— wave channels has been established from 
asymptotic weak coupling analysis in the limit U /W <C 1 
in two and three dimensiona^Ti^j from the observation 
of a positive pair-binding energy computed from ex- 
act diagonalization on various "Hubbard molecules,"— 
from dynamical cluster approximation (DCA)21 and clus- 
ter dynamical mean-field theory (DMFT) calculations^^, 
and from extensive density matrix renormalization group 
(DMRG) studiesi^ii^ of various ladder systems extrapo- 
lated to the thermodynamic limit. 

While there is still some controversy over the strength 
of the pairing tendencies under particular circumstances, 
there is a growing consensus that superconductivity in 
the repulsive Hubbard model is generic under a wide 
range of circumstances. Various physically motivated 
approximate calculations, including numerically imple- 
mented functional renormalization group (FRG)^, dy- 
namical cluster approximation (DCA)^ and cluster dy- 
namical mean- field theory (DMFT) calculations^^, fluc- 
tuation exchange approximation (FLEX)^, Eliashberg^^ 
and self-consistent two-particle calculations^^, as well as 
strong coupling approaches based on variational wave- 
functiona^i and slave-particle mean-field theorie a^^'^° , 
have also strongly indicated that such unconventional 
pairing is present, especially near half-filling, and that 
Tc is maximized in the physically relevant range of inter- 
mediate coupling, U ~ W, and decreases at both larger 
and smaller U. The decrease of Tc at smaller U is due 
to the fact that the strength of any induced attractive 



interaction must vanish as [/ — >■ 0, while the decrease 
as [/ — > oo is due to Mott physics which tends to local- 
ize fermions near particular lattice site a^^i'^° , potentially 
suppressing both the superfluid stiffness^ and the pair- 
ing tendencies^ of the electrons. 

It is broadly (although not universally^) accepted 
that the basic features of superconductivity arising from 
short-range repulsion between electrons are moderately 
generic and must be in play in a broad class of uncon- 
ventional superconductors such as the cupratcs, heavy 
fermion and organic superconductors, Sr2Ru04, and the 
recently discovered Fe-pnictides, despite differences in 
band-structure and local quantum chcmistryj^ In par- 
ticular, estimates of the optimal Tc for d-wave pairing 
in 2D Hubbard models, obtained using a variety of ap- 
proximate computational approache o^^'^'^i^^'^^i^^'^'^i^^'^^ , 
suggest that Max[Tc] ^ 10~^VF/ao where vp is the Fermi 
velocity averaged over the Fermi surface (FS), and oq is 
the interatomic distance. Using vp/a ~ leVr^ one ob- 
tains an estimate of the optimal Tc of order lOOK, com- 
parable to the d-wavc transition temperatures found in 
optimally doped cupratcs. 

However, in addition to the vexing problem of how to 
unambiguously establish or falsify the applicability of a 
particular electronic pairing mechanism to real materi- 
als, at least one key theoretical question remains to be 
addressed: It is well known that longer-ranged compo- 
nents of the electron-electron interaction, even simply a 
nearest-neighbor repulsion V between electrons, suppress 
the pairing tendencies of the Hubbard model. As we will 
discuss below, this can be seen clearly from the struc- 
ture of the asymptotic weak-coupling approach and also 
in exact diagonalization studies of Hubbard molecules. 
This effect has also been investigated in DMRG studies 
of t-J ladders^ (which we extend to Hubbard ladders 
in the present paper). Thus, the issue to be addressed 
is whether the deleterious effects of longer-range com- 
ponents of the electron-electron repulsion make a purely 
electronic mechanism of superconducting pairing physi- 
cally implausible. 



The physics behind the suppression of Tc is transpar- 
ent: an effective attraction which can give rise to un- 
conventional superconductivity in p— wave, d— wave, ex- 
tended s— wave, and other channels emerges in the theory 
from the renormalization (screening) of the local repul- 
sive interaction U by particle-hole fluctuations of the con- 
tinuum of fermions. For small J7, the induced interactions 
are necessarily of order U"^ /W or smaller, and more gen- 
erally one expects the induced interactions to be overall 
weaker than the bare interactions that give rise to them. 
However, if the bare interaction is short-ranged, then it 
only contributes in the trivial s-wave channel, while be- 
cause the polarization of the particle-hole continuum de- 
pends on the energy and momentum transfer between 
initial and scattered fermions, the induced interactions 
are non-local and may have attractive contributions in 
other pairing channels. 

If, however, the bare interaction is extended, so it con- 
tributes a repulsive piece in the same extended s— wave, 
p— wave, d— wave, etc, channels, the competition between 
the bare and induced interactions becomes more seri- 
ous. This issue has been recently addressed in a par- 
ticularly pointed fashion by Alexandrov and Kabanov 
(AK) (Refi^). They have concluded that in models with 
a screened Coulomb interaction, the bare repulsion in 
non-s-wave channels is so significant that it effectively 
overwhelms the induced attractions, thus eliminating as 
physically realistic the entire class of theories in which 
the pairing comes from electron-electron repulsion. 

We will show that the conclusion of Ref<^ is not war- 
ranted. From exact diagonalization and DMRG studies 
discussed below, it is clear that there is no anomalous 
sensitivity of the pairing strength to longer-range inter- 
actions; in an extended Hubbard model, significant sup- 
pression of unconventional pairing occurs only when the 
farther neighbor interactions, V , are a noticeable fraction 
of [/, and not when they exceed the (much smaller) gap 
scale associated with pairing. From the various approx- 
imate extensions of the weak-coupling analysis to inter- 
mediate coupling, it is clear that a key physical ingredient 
missed in the discussion of AK is the near resonant en- 
hancement of the induced attractions at special favored 
wave-vectors when the system is in (not too close) prox- 
imity of a density wave ordering transition; under these 
circumstance, if the favored vectors appropriately nest 
the FS, the appropriate weighted average of the bare and 
induced interactions can favor unconventional supercon- 
ductivity, despite the fact that the induced interactions 
are generically weaker than the bare interactions 

Even in the strictly weak coupling limit, where con- 
trolled calculations can be carried through (as shown be- 
low), the effect of longer range interactions turns out 
to be surprisingly muted. Longer-ranged interactions 
in lattice models (interactions between nearest neigh- 
bors, second-nearest neighbors, etc) generally have com- 
ponents only in some particular pairing channels and do 
not contribute to other channels. On the other hand, 
the particlc-holc polarization bubbles which enter the in- 



duced interactions have components in all channels, and 
quite often more than one component is attractive. In 
this situation, it is quite possible that longer-range inter- 
action either does not affect the leading pairing compo- 
nent at all, or it suppresses the leading attractive com- 
ponent but leaves the subleading one unaffected. 

In this context, we revisit in this paper the issue of 
unconventional superconductivity from electron-electron 
repulsion in an extended Hubbard model on a square lat- 
tice with an onsite repulsion of magnitude U and nearest- 
neighbor and next-nearest-neighbor couplings of magni- 
tudes V and F', respectively. We study this problem 
using the same asymptotic weak-coupling methods that 
were previously applied to the pure Hubbard model, com- 
plemented by a new DMRG study of a two-leg ladder. 
We consider arbitrary fermionic density and analyze the 
pairing problem for various relations between V,V' , U, 
and the bandwidth W. It is important to stress that 
the physics governing U and V is generally quite dis- 
tinct, and they should be considered as essentially inde- 
pendent parameters. For instance, in a transition metal 
oxide, U reflects the atomic physics at short distances, 
while V reflects the screening effects of all the degrees of 
freedom that are integrated out at high energies, includ- 
ing the polarization of the surrounding medium in the 
solid state environment. We find that the interplay be- 
tween the short and longer-range repulsive interactions 
results in a fascinating, complex phase diagram with a 
large number of distinct unconventional superconducting 
phases. The longer-range interactions do indeed tend to 
suppress unconventional superconductivity in some chan- 
nels, but still we find that unconventional superconduc- 
tivity emerges for all densities and for all relations be- 
tween parameters of the model. 

This paper is organized as follows. In the next section, 
we review the weak coupling approach to pairing in a 
system of fermions with generic finite-range interactions. 
In Sec. mil we discuss the case of short-range interactions 
and small electron density, where the fermionic dispersion 
can be approximated by a parabola. In II VI we extend the 
approach to arbitrary electron densities and analyze the 
extended Hubbard model with nearest, and next-nearest 
neighbor interactions. In Sec. |V]we consider the case of 
screened Coulomb interactions, which we can only treat 
approximately, even in the small limit. In section IVTl 
we present the results of the DMRG calculations of the 
pairing amplitudes in Hubbard ladders. We discuss the 
results and present our conclusions in Sec lVIII 



II. ASYMPTOTICALLY EXACT 
WEAK-COUPLING APPROACH 

The basic idea of the Kohn-Luttinger (KL) approach to 
superconductivity is that the pairing interaction is given 
by the irreducible vertex in the particle-particle channel, 
which is generally different from the bare interaction be- 
tween the two given fermions and includes contributions 
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from the continuum of particle-hole excitations. These 
additional contributions to the pairing interaction can 
give rise to an attraction in a particular pairing channel 
even if the original interaction between two fermions is 
uniformly repulsive. 

The KL approach can be formally justified if there 
is a separation of scales: superconductivity is assumed 
to come from fermionic states very near the FS, while 
the renormalization of the irreducible pairing interac- 
tion comes from fermions with energies comparable to 
the Fermi energy. In this situation the fully renormal- 
ized pairing interaction can be approximated by its value 
when both incoming and outgoing fermions are on the FS 
and the two incoming fermions have opposite momenta 
and zero energy (where the energy is measured as devi- 
ation from the chemical potential). Such a separation 
of scales clearly occurs when the interaction U is much 
smaller than the fermionic bandwidth W; under appro- 
priate circumstances, it may hold approximately (in a 
physical sense) even when U is comparable to the band- 
width2^. 

The KL approach can be implemented in terms of a 
two-step renormalization group (RG) procedure^: In the 
first step, we integrate out all modes outside a narrow 
range of energies Oq about the Fermi energy. Hq is not a 
physical energy in the problem, but rather a calculational 
device. It is chosen large enough so that the interac- 
tions can be treated perturbatively (e.g., the Cooper log- 
arithms arc still small), but small enough that it can be 
set to zero in all renormalizations in non-singular chan- 
nels without causing significant error, i.e. it is chosen to 
satisfy the inequalities 



NpU^ > J7o > EFexp{-[l/{NFU)]}, 



(1) 



where Np is the density of states (DOS) at the Fermi 
energy (for both spin orientations) and Ep is the Fermi 
energy. The effective interactions generated in the pro- 
cess then serve as input interactions in a second step, in 
which the remaining problem is solved using the logarith- 
mical RG techniqu o^^i^'' which is equivalent to BCS when 
only the pairing channels are singular. (It involves a more 
complex parquet RG if some other channel, e.g., SDW or 
CDW, is also singular and competes with superconduc- 
tivityjiii^l a situation we will not consider explicitly for 
the present.) Tc is, up to a multiplicative constant, given 
by the energy scale, T*, at which the pairing interaction 
grows to be of order one. It was shown by explicit pertur- 
bative calculation of Tc up to 4th order in U that the re- 
sulting expression for T* is independent of fio (Refi^i^). 
The generic prescription for computing the leading or- 
der asymptotic behavior of Tc for weak interactions is 
the following^.: first, compute the effective interaction in 
the Cooper channel at energy scale Qq perturbatively, to 
second order in the bare interaction. The corresponding 
diagrams are shown in Fig. [TJ We denote this effective in- 
teraction as r'"' (k, — k; p, — p) = F^"' (k, p), where k and 
p denote points on the FS and the superscript a differ- 
entiates between whether the electron pair forms a spin 



k, (7 



fc'.cr 



(2a) 



(2b) 



(2c) 



(2d) 



(2e) 



FIG. 1: First and second order diagrams which contribute 
to tlie effective interaction r"(k, p) in the Cooper chan- 
nel. The solid line corresponds to an electron propagator 
whereas the dashed line represents an interaction vertex, a 
and a' axe spin indices. The singlet and triplet components 
r'^' and r'*' are obtained by antisymmetrizing the interac- 
tion and taking symmetric and antisymmetric combinations, 
(r(k, p) + r"(k, -p))/2 and (r(k, p) - r''(k, -p))/2, respec- 
tively. 



singlet (r*^**^) or a spin triplet (F*^*^). Wc then construct 
the related dimensionless matrix 



7('') (k, p) = Nf v/5F/t.F(k)F(") (k, p) ^vf/vf{v), (2) 

where i'F(k) is the magnitude of the Fermi velocity at a 
given point k on the FS, and vp and Np are the average 
Fermi velocity and average DOS at the FS. For a circular 
FS, vp{)<i) = vp = vp and Np = Np, but the approach 
also holds for lattice models of fermionic dispersion. 
Since 7^°^(k, p) is a real, symmetric matrix, it has a 

complete set of eigenstates 4'^''''' and eigenvalues A*^"''^ 



(3) 



Since all the systems we will consider are inversion sym- 
metric, the spin symmetry is implicitly determined by 
whether / transforms under an irreducible representation 
that is even or odd under inversion, so wc will henceforth 
leave it implicit, (a, I) {I). Among all the possible 
solutions, we identify the most negative eigenvalue. 



A = Min 



A < 



Then, 



T,^^i.exp[-I/|A|]. 



(4) 



(5) 



Since the effective interaction F is computed perturba- 
tively assuming weak interactions, A can be expressed as 
a power series in the bare interaction. Assuming that the 
interaction is a local Hubbard U, wc have 



A 



[l + ai|u| + a2i 



(6) 
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for negative U {u = NpU), and 

A = -M^ [60 + biu + b2i 



(7) 



for positive U. ai and bi are non-universal constants 
which depend on diniensionahty, crystahine point group, 
electron concentration, and details of the banstructure. 
The absence of a linear in u term in A for [/ > is 
a consequence of the fact that in a repulsive Hubbard 
model, an attractive interaction in any pairing channel 
can only appear due to a renormalization of the original 
interaction by the particle-hole continuum. As a conse- 
quence, Tc is a highly non-analytic function of u in the 
limit u — > 0, with different functional dependencies for 
positive and negative u. 

In the next two sections we apply this general proce- 
dure to systems with isotropic dispersion, and to systems 
with full lattice dispersion. 



III. CONTINUUM ELECTRONS WITH WEAK 
SHORT-RANGED INTERACTIONS 

Dilute electrons on a lattice treated in the effective 
mass approximation are equivalent to the continuum 
problem in which the FS is a sphere (a circle in 2D) 
and the electronic dispersion has the parabolic form 
Efc = k^/2m — fi. As a first example, we thus con- 
sider the low density limit of electrons with an inter- 
action, U(r) of finite range, tq, which is independent 
of the electron density. There are thus two indepen- 
dent small parameters: kpTQ <ti 1 and kpa <^ 1 where 
a = m[/(0)/47r ^ mU{r — 0)rg is the s-wave scattering 
length in the Born approximation, and U is the Fourier 
transform of U . 

To first order in U, the singlet and triplet components 
of the pairing interaction in the momentum space are 

r(«) (k, p) = [;7(k - p) + ;7(k + p)] /2 
r(*)(k,p)= [[/(k-p)-c7(k + p)] /2 (8) 

For k and p located at the FS, (k - p)^ = 2fc|,(l-cos6'). 
Since ?7(k) changes slowly over the range of momenta 
< |k| < 2fci?, it can be expanded in powers of kpro <C 1: 



U{k) = U{0) (1 + ai\k\^rl + a2\k 



(9) 



where a„ are pure numbers related to moments of U{r). 
The first term in ([9]) then contributes to s— wave pair- 
ing, the second term contributes to s— wave and p— wave 
pairing, the third term contributes to s— wave, p— wave, 
and d— wave pairing, and so on. Indeed, because there is 
full rotational symmetry, the pairing wave functions (i.e. 
the eigenstates in Eq. [3]) are completely determined by 
symmetry [i.e. they are uniquely specified by angular 
momentum sectors labeled by /). Correspondingly, the 
eigenvalues, A^'-*, can be computed directly as the appro- 
priate symmetry determined average of 7 over the FS. In 



3D, Np = mkp / {2tt'^), and it is important to note that 
u = NpU{Q) = 2{akp)/TT . 

Then, following the renormalization group procedure 
outlined above, explicit perturbative expressions can be 
obtained for A^'^ which, to second order in m, are 



A^'' = {kpa)A^^'^ + {kpafB^^^ + , 



(10) 



where A^'^ and i?''^ are dimensionless functions of the 
dimensionless variables, kpr^ and Npflo (where f^o is 
the scale at which the couplings are defined, as discussed 
above Eq. [T]). However, we are interested only in their 
leading behavior for small values of both these parame- 
ters. Starting with the first order terms, A^'^ is mani- 
festly independent of Npflo but varies as kpro — >■ ac- 
cording to A^') ^ [kprof^. Specifically, for the three 
lowest angular momentum channels. 



AW = ^{rokpf . 



[rokpf, 



(11) 



The components with larger I are progressively smaller 
in powers of r^kp. To this order, the sign of A^'-* depends 
on the sign of a and on the sign of a/. For a repulsive 
Yukawa potential (i.e. for tj{q) = {/(0)/[l + [qrof]), 
a > 0, q;i < and a2 > 0, and this is rather generic for 
purely repulsive interactions. Thus, there is no Cooper 
instability to first order in [kpo). 

The situation is altered by corrections to r(")(k, p) to 
second order vaakp, see Fig. [TJ The diagram 2a describes 
the renormalization in the particle-particle channel which 
makes a contribution to i?*-*' in the s-wavc channel pro- 
portional to log[Vl^/r2o]- However, so long as the inequal- 
ity in Eq. [1] is satisfied, this only produces a small cor- 
rection to A*-*^ relative to the first order repulsive term. 
The other four second-order diagrams which contribute 
to i?*-'^ contain particle-hole bubbles and account for the 
important renormalizations of the irreducible pairing in- 
teraction. 

A key result obtained already by KL (Refill) is that, 
as kpTQ — > 0, -B*^'' Pi ^ Q. KL further showed that 
for sufficiently large I > 0, S*-'^ are negative, and that 
while A^''' falls exponentially with increasing I, B*^') ^ 
l^^. From this, they concluded that superconductivity 
in some (possibly high I) channel was inevitable. 

We focus attention on the physically more important 
cases of relatively small I. It has been shown— 1^ that /?/ 
are negative for ^ > 0. As a result, 

ttA'p' = {akp) [\ai\{rokFf ~ \Pi\{akF) + Oikpaf] , 
TrAf"^) = [akp] [\a2\{rokpf - \l32\{akp) + ...], (12) 

where all /3; are of order 1. The issue then is which 
terms are larger, the negative second-order or the pos- 
itive first-order contributions. Clearly for kp^ ^ ^ a, 
the repulsive part of the interaction is dominant. How- 
ever, there is a broad range of circumstances in which 
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In this situation, the attractive interac- VW/U'^ = a > held constant), and V = a'U (i.e. 



tions dominate, even for relatively small / > 0. In 3D 



/3i = - 



8(21og2- 1) 
5^ 



P2 - 10-2/3i 



(13) 



(see Refs. 013) > and hence at small density a 3D Fermi 
system with short-range repulsive interactions undergoes 
a wave pairing transition with 



TcOcexp[-l/{\Pi\{akFf)] 



(14) 



(For the full calculation of Tf, sec Rcf4i.) In 2D, the cal- 
culations arc a bit more tricky because the static particle- 
hole bubble for free fermions is independent of momen- 
tum transfer for|k — p| < 2k p, so one has to go to the 
next, third order, to obtain the momentum dependence 
of the dressed interaction (seeRefQ). The final result is, 
however, similar to the one in 3D: a 2D Fermi system at 
a small density undergoes a p— wave superconductivity 
with a low Tc- 

Note that a similar analysis was carried out in Ref |40| 
for the extended Hubbard model at small electron densi- 
ties, where instead of a Yukawa interaction, strong Hub- 
bard interactions were considered with U ^ V ^ the 
bandwidth with similar results to those discussed here. 



IV. EXTENDED HUBBARD MODEL AT WEAK 
COUPLING 

As a second more directly experimentally relevant ex- 
ample, we apply the reasoning from Sec|TT] to a higher 
density of electrons in the extended Hubbard model on 
a two dimensional square lattice with the Hamiltonian: 



H 
Ho 



Hq + Hint 



(15) 



Hi, 



U ni^riii + V TiiTij + V' niUj 



(16) 



where Cio- is the destruction operator of an electron on lat- 
tice site i, with spin a, and ((«, j)) denotes nearest- 
neighbor and second-neighbor pairs of sites respectively, 
rii = c\^Cia and the average density of electrons per 
site is ri =< Ui >. For simplicity, we consider band 
electrons with only nearest-neighbor hopping. As is well- 
known, such a system possesses a non-generic particle- 
hole symmetry. This feature, however, does not play a 
significant role in the resulting phase diagram since we 
consider the model away from half-filling. (See, however, 
Ref. [13). 

To begin with, for = we show that the phase 
diagram is drastically different depending on the ratio 
F/C/, even as u ~ U /t 0. Different asymptotic analysis 
is required for the case V — aU'^ /W {i.e. for w — !• with 



for u with V/U = a' > held fixed). We find 
that in both cases, the ground state is superconducting 
for all dopings, but the symmetry of the pairing state is 
generally different. We then add V ^ V and show that 
it adds other pairing states to the phase diagram. 



A. V V^/W, y = 

Since the effective interaction is computed to 0{U'^), 
when V = aU"^ /W , we only need take V into account to 
first order (Diagram 1 in Fig. [1]), i.e. at the bare level. 
A non-zero V then produces a correction to the effective 
r(k, p) in the form 

'5r(k,p) = a— [cos {hj; - pr,) + cos [ky - py)] 



"l7E^^''''nk)^("-^)(p) 



W 



(17) 



where in the second line the sum is taken over appro- 
priate basis functions defined on nearest-neighbor sites 
with Aig or extended s-wave symmetry, Big or da;2_,y2- 
wave symmetry, and E„ or p-wave symmetry: 



Iff 



A 

B2g ■■ 

Eu '■ 



0(^-1) (k) = [cos(fc,) +cos(fcy)]/V2, (18) 
0(^'-y'.i)(k) ^ [cos(fc.O - cos{ky)]/V2, 
(^("■i)(k) = sin(fc,), 0(^^i)(k) = sin(fcj,). 



Thus, the nearest neighbor interaction acts as a sepa- 
rable repulsive interaction within these subspaces. Be- 
cause all three components are repulsive, V tends to 
suppress the pairing tendency in all three of these chan- 
nels. However, note that unlike the case in the contin- 
uum, there are multiple (infinite) orthogonal functions 
defined on the FS which transform according to each ir- 
reducible representation of the point group - for instance 
0(s^2) ^ [cos(fc^) + cos(fcj,)]^ also transforms according to 
Aig under operations of the point group. Thus, even for 
large a, the presence of a repulsive first order term in a 
given channel does not preclude the existence of a more 
"extended" form of pairing in the same channel. 

It has been shown previously in various studies of the 
repulsive U Hubbard model that near half-filling, ^n^n^idl 
predominant pairing instability to order C/^ is in the 
da;2_j,2 channel, and the subdominant pairing eigenvalue 
occurs in the or g-wave channel, while somewhat 
further from half-filling, at n < 0.62, there is a range of 
electron concentrations for which the B2g or dxy pairing 
solution is dominant: Representative gap functions (of 
the shortest possible spatial range) with these symme- 
tries are: 

0(9.3) ^ [cos{krc) - cos{ky)] sin{kx) sin{ky). (19) 
0(^2/.2) _ sin(fc^)sin(fcj,). 
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FIG. 2: Phase diagram as a function of electron concentra- 
tion for the 2D extended Hubbard model in the regime where 
V ^U'^jW and V = 0. The region of extended s— wave su- 
perconductivity is labeled as 1. For F = 0, our phase diagram 
is similar but not identical to the one obtained by Hlubinai^. 
In his calculation, there is a small region of p— wave super- 
conductivity between dxy and d,j.i_yi phases. We found only 
d-wave states. 
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FIG. 3: The effective interaction, Ve// = X/Nf where Np 
is the density of states at the Fermi level, as a function of 



V'/V = Q2. 

non-zero V' . 



(n = 0.6) with V = aUyW , a = 0.16 and 



Figure [2] shows the phase diagram to order , which 
wc obtained numerically for V = all'^/W, as a func- 
tion of a and n. Since the analysis considered here 
breaks down for concentrations sufficiently close to half- 
filling because of competition with antiferromagnetism, 
we have investigated solutions only for < n < 0.95. 
For 0.76 < n < 0.95, a finite a > 0.2 destabilizes the 
dx'i-y'i solution in favor of the subdominant g-wave solu- 
tion (labeled as xy{x'^ ~ J/^) on the phase diagram in Fig. 

For 0.68 < n < 0.76 a finite a again destabilizes the 
dx2-y2 solution, but the state that emerges instead has 
a particular extended s— wave symmetry, that is domi- 



nantly of the form Re(a; -|- To shorten notations, 

have labeled it as "1" . This state has zero amplitude for 
on-site pairing in this limit but transforms nevertheless as 
a trivial irreducible representation of the tetragonal point 
group. Note that since the wave state is Im(a:: + iy)'^, 
in the continuum limit, the "1" and g states are degen- 
erate corresponding to angular momentum £ = 4 pair- 
ing. However, lattice effects lift this degeneracy. Indeed, 
there is a phase transition between wave and extended 
s— wave state at n = 0.76, which reflects a level crossing 
of these two eigenvalues, both of which are sub-leading 
at y = 0. Note that this phase boundary is vertical since 
the two solutions are unaffected by V to first order. For 
0.58 < n < 0.68, we have found that a finite a favors 
the dxy solution. The phase boundary between d^y and 
extended s— wave phases is also vertical. 



B. 



Next, we consider the effect of a non-zero second neigh- 
bor interaction V' on the phase diagram. Specifically, we 
take V = a2V,0 < < 1. In this case, V - U'^/W, 
so like V, its effects can be computed via the first order 
diagram 1 in Fig. [TJ The expansion of the V' interaction 
into angular harmonics is straightforward, and in addi- 
tion to terms already present in Eq. (|17p . one finds the 
following contribution to the effective Cooper channel in- 
teraction from V': 



Srv'{k,p) = 2aa2—cos{kx -px) cos {ky - py) (20) 



This term is a sum of separable repulsive interactions in 
the extended-S (Aig), dxy (B2g), and p-wave (E„) sub- 
spaces. It is important to stress that it does not affect the 
da;2_j^2 (Big) or the g-wave (A2g) subspaces. Thus, rel- 
atively close to half-filling where the leading eigenvalues 
belong to the dx2-y2 and g-wave subspaces, the phase 
boundaries are unaffected by V . We have studied the 
phase diagram in the presence of non-zero V' and have 
found that while the phase boundaries are affected quan- 
titatively by non-zero a2, the topology of the phase di- 
agram itself is unchanged. We have observed that the 
pairing strengths are surprisingly robust in this regime, 
as a2 is increased. Figure El displays the effective interac- 
tions Veff = \\\/Nf, where Np is the density of states at 
the Fermi energy, in the presence of non-zero V' . For dcf- 
initeness, we work in a fixed density of electrons n = 0.7, 
and fixed a = 0.16 and show the effective interaction 
as a function of a2- The dx2-y2 pairing strength is the 
weakest one here and is not shown. It is apparent from 
the figure that while the pairing strengths of the dxy and 
extended s-wave states are affected by non-zero a2, this 
repulsion is not strong enough to suppress superconduc- 
tivity altogether. 
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C. V r^U 

Next, we consider the case when V ^ U, in which case 
the bare repulsive interactions are of the form 



r7(fc,fe') = U + Vg{k- k') 



(21) 



where the subscript on F above is to remind the reader 
that this comes from diagram I in Fi^ The dimension- 
less function 

g{q) = cos qx + cos qy + 2f3 cos q^ cos qy (22) 

specifies the momentum dependence of the interactions, 
where /3 = V'/V. We must now compute the perturba- 
tion expansion of the effective interaction in the Cooper 
channel to second order in all the bare interactions which 
has the form 

F(^)(fc, k') = Ti{k, k') + U^f['\k, k') + V^fi'\k, k') 
+UVfi'\k,k') + --- (23) 

where s denotes the spin-singlet channel. It is straight- 
forward to show that 



f['\k,k') = Npln 



W 



X{k + k') + 0{no) (24) 



where Np is the density of states at the Fermi energy. 
The first term is obtained from diagram 2a and the sec- 
ond from 26 in Fig. [TJ Similarly, 



f^''\k,k')=NF{g*g) (fc-fc')ln 



W 



Xi{k,k') 



g{k - fc')X2(fe, k') - 2 [g{k - fe')] X(fe - k') + O(17o) 

(25) 



where the first term comes from diagram 2a, the sec- 
ond from 26, the third from 2c, 2d and the last term 
from 2e in Fig. [TJ The convolution that enters in the 
first term is found to be (assuming that the FS averages 
< cos 2/3; >=< cos2/.y >= 0, where 1 and —1 are inter- 
mediate fermionic momenta in Fig. [1^) 



- (cos qx + cos qy) -|- cos q^ cos qy. (26) 



Lastly, 



/^"^(fc,fc') = aNFg(k-k')\v, 



W_ 



-X2{-k,k') 

+ X2{k, k') + 2g{k - k')x{k - k') + 0{VL^) 

(27) 



where a is a constant of order unity. The first term is 
obtained from diagram 2a, the second from 26, the third 
from 2c, and the last from 2e in Fig. [TJ The functions 
x(q), Xi(fe, k'), and X2{k, k') are generahzed susceptibil- 
ities and are all expressible in terms of the one electron 
Matsubara Green function, G{p) = {iujp — Cp)~^'- 



G{p)G{p + q), 



{2n)'i+^ 

Xi{k,k') = [ g{k-p)g{p-k')G{p)G{p-k-k') 
Jp 

X2{k,k') = f [g{p + k)+g{p^k')]G{p)G{p + k-k') 
Jp 



(28) 



Similarly, in the spin-triplet channel, the effective inter- where 
action takes the form 



rW(fc,fc')= Vg{k-k') + U'f['\k,k') + V'f^'\k,k') ^^^^ 



2 At), 



UVf^'\k,k') 



fi'\k,k') = -xik-k') 
A'\k,k') = fi'\k,k') 
f^\k,k') = -2g{k-k')x{k-k') 



(30) 



P ~ the contribution from diagram 2a disfavors 
(29) solutions which involve nearest-neighbor pairing. For a 
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risks double-counting of certain important physical pro- 
cesses, which are already implicit in the screened form 
of the interaction. Thus, to treat the Coulomb problem, 
we are forced to adopt a more qualitative, although still 
physically sensible approach, in lieu of the asymptotically 
exact approach we have pursued up until this section. 
For simplicity, we again consider the explicit case of di- 
lute electrons, in which the dispersion is quadratic. It 
is important to note, however, that for small r^, where 
weak-coupling intuition is most likely to be correct, the 
screening length is parametrically larger than kp^, so 
even the screened interaction is not, in any naive sense, 
short-rangedjSii^ 

In 3D, the screened Coulomb interaction within the 
RPA approximation is 



FIG. 4: The effective interaction K// = X/Nf for the g-wave 
state. Here, we have chosen V — 0.5U and V' — 0.5V . All 
other eigenvalues are subdominant and are not shown here. 



non-zero /3, solutions with second-neighbor pairing am- 
plitude are also suppressed. 

In Fig EJ we show the pairing strength of the A2g (g- 
wave) state when V = 0.5U, V = 0.5V. We have found 
that the pairing strengths in other channels are signifi- 
cantly weaker and are therefore not shown in the figure. 
We note that, interestingly, the g-wave pairing strength 
is much greater in the regime V ^ U since it gains en- 
ergy from the attractive effective interactions that are 
obtained at 0{V'^), which were neglected in the regime 
where V ^ U^. Thus, while the states with nearest 
neighbor pairing are severely suppressed by momentum- 
dependence of the first-order terms, states involving fur- 
ther ranged pairing are more enhanced by second-order, 
KL contributions than they are suppressed by first-order 
contributions. It is important to stress that in the weak 
coupling limit, such enhancements cannot be ascribed 
to a well-defined bosonic 'glue', such as those associated 
with proximate spin or charge density states, since insta- 
bilitites toward spin and charge ordered phases occur at 
weak coupling only in exponentially narrow regions in pa- 
rameter space (close to a perfectly nested Fermi surface, 
for instance). 



V. CONTINUUM ELECTRONS WITH 
SCREENED COULOMB INTERACTIONS 

The Coulomb interaction is long-ranged, and so never 
can be safely treated in perturbation theory. Physically, 
in a good metal, it should be screened, and so equivalent 
to an appropriate finite-range model, but because screen- 
ing at long distances involves electrons with energies ar- 
bitrarily close to the Fermi energy, this is something that 
is subtle to treat in a controlled RG. Simply to replace 
the Coulomb interaction with a screened Coulomb inter- 
action, and then to carry out the usual RG from there 



Uiq) 



47re^ 



where 



K = O.Slkprl^^ , and = 1.92e^ /hvp. If we calculate the 
s-wave scattering length and the range of the interaction, 
as wc did in Section lllll then akp = rne^kp/ ~ 0.8 and 
To ~ 1/k ^ a/^rs Y^^, so in the "weakly interacting limit" 
rs << 1, it follows that tq ^ a. Because akp — 0{1), 
an expansion in powers of the interaction is not possi- 
ble. In addition, the screening of the Coulomb interac- 
tion already contains one of the diagrams (a particle-hole 
bubble) which for short range interactions participated in 
the KL rcnormalization of the irreducible pairing inter- 
action. None-thc-lcss, it is still plausible that the same 
sort of KL calculation (with the screening diagram omit- 
ted) will still give a physically reasonable account of the 
pairing in this case, although the formal justification for 
this approach is on somewhat less secure footing. 



4^1 



Skpq 



q- 2kF + q 
log- 



2kF-q 



(31) 



(32) 



A. Ts « 1 

For rs << 1, since the relevant values of the mo- 
mentum transfer which enter the various calculations is 
|q| ^ K « kp, lJ{q) can be well approximated by taking 
n(g) w n(0) = 1. Then 



U{q) = 

and consequently 

NptJiq) = 0.041 



47re 



1- cose + 0.lMrs 



(33) 



(34) 



where q = ^/2kp[l — cos(6')]. Expanding this bare U{q) 
in angular harmonics one finds that they are all posi- 
tive (repulsive), but that their magnitudes decay with 
increasing /, eventually as e 



-Irt 



for Ir^^ » 1 (Ref, 



l9^ 
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The KL-type contributions to A*^'-* come from particle- 
hole renormalizations which involve particlc-holc bubbles 
(diagrams 2b, 2c and 2d in Fig[T]). For ^ > 1, typ- 
ical momenta associated with these bubbles are of or- 
der kp, and since NpU{kp) = 0{rs) ^ 1, this implies 
that these terms make contributions to A*-'^ which are 
small compared to the scale of the repulsive contribu- 
tions ^ NpU{0) = 0(1). Obviously then, for I not too 
large, A^'^ > 0, i.e., there is no superconductivity. The 

1 /2 

situation is different for large / > l/rs owing to the fact 
that the exchange diagram (diagram 2b in Fig|T]) has a 
non-analytic dependence on 9 near 9 = 0, i.e. when k 
and p are nearly parallel. Moreover, under these circum- 
stances, the typical internal momenta in the particle- hole 
ladder are also nearly parallel to k and p, from which it 
follows that the relevant interactions involve near-zero 
momentum transfer where the interactions are strong, 
NfU{0) = 0(1) (Ref.iS). Because of this smaU 9 non- 
analyticity, the contribution to A*-'-* from the exchange di- 
agram scales as 1/Z^ rather than exponentially in I, and so 
makes the dominant contribution to A*-'-* for large enough 
I. This non-analytic contribution is attractive indepen- 
dent on the parity of I (Refii). Of course, Tc is always 
very small for such high I pairing, but, these arguments 
suggest that KL-type superconductivity survives. 



B. r^c > Vs > 1 

For Ts ^ 1, even in the continuum, electrons form a 
Wigner crystal state, which is manifestly neither a Fermi 
liquid nor a superconductor. However, since the critical 
value, r^ijc, for Wigner crystallization appears to be nu- 
merically large, it is reasonable to wonder what happens 
when Ts is in the range 1 < < r^c- In this range, there 
is no guarantee that the Fermi liquid fixed point is even a 
correct starting point for examining the problem - there 
could be a variety of possible other phases, and, even if 
one is prepared to ignore all standards of mathematical 
rigor, there is no controllable expansion because NpU{q) 
is of order one for all q which connect points on the FS. 

The approach adopted by AK to address this point 
was to assume the same perturbative approach that we 
have outlined for <C 1 can be applied for > 1, i.e. 
the pairing vertex is calculated to second order in the 
screened Coulomb interaction which is taken to have the 
Yukawa form as in Eq l33l and only the exchange diagram 
2e from Fig. [T]is evaluated to determine the second-order 
contribution to the pairing vertex. Evaluating A*-'-* for the 
first few angular harmonics, one then finds^ that the in- 
teractions in the p and d— wave channels (among others) 
are repulsive for Ts 1, and only become attractive again 
for unphysicaly large values, i.e., > 54 for p— wave and 

> 26.5 for d— wave pairing. From this AK concluded 
that the KL mechanism is not a plausible mechanism of 
superconductivity in systems with Coulomb interactions 
with ~ 1 or greater. 

Given that even the starting expression for the 



screened Coulomb interaction in Eq. |3T]makes little sense 
in the large rg limit, where the Thomas Fermi screening 
length is parametrically smaller than the distance be- 
tween electrons, it is difficult to judge the validity of this 
conclusion. We could, with no less justification, repeat 
the AK calculation but keeping the full RPA expression 
for the screened Coulomb interaction, Eq. [211 rather than 
the Yukawa form in Eq. 1331 In this case, an attractive 
pairing vertex is obtained already in first order! (Second 
order contributions from diagrams 2b-2d in Fig. [T]would 
further enhance the pairing tendencies.) As was shown in 
Ref. [H, the largest attractive pairing component attrac- 
tive component is Z = 1 for which A(/ = 1) « —0.07. A 
very similar X{1 = 1) w —0.06 has been obtained in Ref4^ 
where the effective interactions in spin and charge chan- 
nels were chosen to fit Monte-Carlo data for compressibil- 
ity and spin susceptibility. In the absence of a justified 
theoretical computational scheme, it is difficult to gauge 
which approach is better. It is likely that wave is the 
leading attractive pairing component at >> 1 but at 
which Tg it prevails and what is Tc are the two issues not 
settled yet. 

There are other issues in the large limit. Due to 
the proximity of the Wigner crystal, the effective interac- 
tions are strongly momentum-dependent. Magnetically 
ordered states are also natural in this limit, and proxim- 
ity to them can give rise to structure in the pairing vertex 
related to magnetic fluctuations: This observation un- 
derlies theories of p— wave pairing near a fcrromagnetio^ 
and d-wavc pairing near an antiferromagnetic instabil- 
ity^. However, the quasiparticles themselves may be 
largely incoherent, making it difficult to know precisely 
what is being paired. Even if Fermi liquid theory remains 
valid, strong forward scattering interactions, which are 
likely to occur in this limit, can produce large changes 
in the Fermi liquid parameters, with major consequences 
for any theory of superconducting pairing. 



VI. DMRG SOLUTION OF EXTENDED 
HUBBARD LADDERS 

To go beyond the weak coupling limit, where the en- 
ergy scales associated with superconductivity are expo- 
nentially small, one has to resort to numerical methods. 
In order to address the effect of extended interactions 
in the intermediate coupling regime, we have performed 
DMRQi^ simulations of the extended Hubbard model 
(Eq. [ISD on ladders of size L x 2. For V = V = 0, 
this system is knowoi^i^"— to have a spin gap in the 
thermodynamic limit when the electron density is close 
to 71 = 1. The superconducting correlations are d-wave 
like, in the sense that the pair amplitude has an oppo- 
site sign on x and y oriented bonds, and falls off with 
distance as a power law which depends on n and U/t. 
Here, we examine the sensitivity of the superconducting 
tendency of the two-leg ladder to adding extended in- 
teractions. We do this by calculating the spin gap and 
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FIG. 5: The spin gap As extrapolated to the thermodynamic 
limit for two-leg ladders with electron density n = 0.9375 
(circles) and n = 0.875 (diamonds) and t = 1, U = 8 and 
V = 0, as a function of V. Beyond V ~ 2.6, a transition 
to a state with pronounced charge density oscillations at a 
wavevector close to (tt, tt) is observed. 



the decay of a proximity induced superconducting order 
parameter as a function of distance, for various values 
of V, v. This study extends earlier results on the t-J 
model with a nearest-neighbor V;^ which found that the 
pairing survives up to F « 4 J. 

The spin gap is defined as = E{S = 1) - E{S = 0) 
where E{S) is the ground state energy with spin S, ex- 
trapolated to the thermodynamic limit L — >■ cx), as a 
function of V. The system sizes used in the extrapola- 
tions where L = 16, 32, 64. The spin gap is then extrap- 
olated to L — > oo by fitting As(l/L) to a second order 
polynomial. The extrapolated value of is up to a fac- 
tor of 2 smaller than As{L = 64); this extrapolation is 
the largest source of error in our results. The following 
parameters were used in the calculations: t = 1, U ^ 8, 
V' = and n = 0.875,0.9375. As V increases, the spin 
gap decreases gradually from its = value, but re- 
mains finite up to y « 2.5. 

For larger values of V, we observe a transition to a 
charge density wave (CDW) state, in which there are pro- 
nounced oscillations in the electron density at a wavevec- 
tor close to Q = (7r,7r). We have checked that the CDW 
transition occurs at 1^ ss 2.5 even for the undopcd (n = 1) 
system. The doped CDW state supports gapless spin 
excitations at the edges, but has a large bulk spin gap 
As « 0.1. (We infer this by noticing that in the lowest 
triplet excitation with Sz = 1, (Sz) is non-zero only close 
to the edges. We have also computed the gap to an exci- 
tation with Sz = 2, in which (Sz) is non-zero in the bulk, 
and found that this gap is finite in the thermodynamic 
limit.) 

We have also examined the effect of a second-neighbor 
V' on the spin gap. Fig. [5] shows the spin gap as a 
function of V tor V = 1, n = 0.9375. Again, we find that 



0.035 
0.03 

0.025 
0.02 

0.015 
0.01 



n=0.9375 



0.2 0.4 
V 



FIG. 6: The spin gap As{L — oo) as a function of V' for 
two-leg ladders with n = 0.9375, t = 1, U = 8 and V = 1. 
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FIG. 7: Induced superconducting order parameter as a func- 
tion of position in a calculation with an edge pair field 
(Eq. I36|) of magnitude A = 0.5. In these calculations, 
t = 1,U = 8,71 — 0.9375. Solid (dashed) lines correspond 
to L = 32 (L = 48), respectively. Results for V = 0,1, 2.5, 3 
are shown. The inset shows the induced order parameter on 
nearest and second-nearest neighbor bonds near the middle 
of an 1/ = 48 system with V = 1. 



while As decreases monotonically upon increasing V' , its 
effect is not dramatic. E.g., upon reaching V = 0.5t, Ag 
has decreased to about 50% of its F = 1, y = value. 

The superconducting response of a ladder system can 
be characterized by the rate at which an externally 
induced superconducting order parameter at the edge 
decays as we move into the bulk. In a gapless one- 
dimensional system, this amplitude decays as a power 
law; in a two-leg ladder with a spin gap, this power law 
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can be shown to be equal to ^ where Kc is the Lut- 
tinger parameter of the (gapless) even charge mode. Fig. 
[7] shows the induced order parameter on a y bond, 

Py{A = [cT(a^, l)ci(x, 2) - c;(x, l)ct(x, 2)]), (35) 

as a function of position a;, on a log-log scale. In this 
calculation, the following boundary term was added to 
the Hamiltonian: 

i/edgc = A [ct(l, l)c^(l, 2) - c^(l, l)ct(l, 2) + H.C.] , 

(36) 

with A = 0.25. This term mimics a proximity-induced 
gap at the edge due to a nearby bulk superconductor. 
The results are shown for systems of length L = 32,48, 
density n = 0.9375, and F = 0, 1, 2.5, 3. 

The induced superconducting order parameter is seen 
to decrease monotonically upon increasing V . However, 
the slope of log(|Py|) vs. log(a;) (the power with which 
the superconducting order parameter decays) far away 
from the edge is not strongly dependent on V ^ except for 
V > 2.5. Although our systems are not long enough to 
allow an accurate estimate of the slope, one can roughly 
estimate Kc ^ 0.4—0.6, well within the range of divergent 
superconducting correlations Kc > 0.25, and close to the 
value in which superconducting and CDW correlations 
decay with the same exponent, Kc = 0.5. 

The inset of Fig. [7]shows the induced superconducting 
order parameter on various bonds near the middle of the 
L = 48 system with V = 1. As can be seen in the figure, 
the order parameter is "(i^2_j,2— like", in the sense that 
the pairing amplitude on x and y oriented bonds is op- 
posite in sign, and the amplitude on the diagonal (next- 
nearest neighbor) is relatively small. The order param- 
eter has a (1^2 _y2— like structure for0<l/<V^R:!2.5. 
For larger values of V (in the CDW phase), the order 
parameter has extended s-wave structure, in which the 
pairing amplitude has the same phase on x and y ori- 
ented bonds. In this regime, however, the order parame- 
ter is much smaller than for V < Vc and it decays much 
faster as a function of distance from the edge (possibly 
exponentially) . 

VII. DISCUSSION 

In this paper we have shown, using a variety of meth- 
ods, that unconventional superconductivity arising di- 
rectly from electron-electron interactions can survive 
even in the presence of longer (but finite) ranged inter- 
actions. Firstly, we have shown that the Kohn-Luttinger 



effect survives even when the range of the interaction 
exceeds the s-wave scattering length, provided that the 
range ro satisfies ro <^ kp^. In this regime, the leading in- 
stability need not occur for an unphysically high angular 
momentum and in most cases occurs for i = 1. For lat- 
tice electrons near half-filling in the weak-coupling limit, 
we have found that the dx2_y2 superconductivity that 
occurs in the Hubbard model survives in the presence of 
longer ranged repulsive interactions V^V ■ ■ , provided 
that these interactions are no larger than aU"^ /W , where 
a is a constant of order unity. Lastly, we have shown, 
using DMRG calculations, that in the intermediate cou- 
pling regime, the spin-gap survives even in the presence 
of a substantial nearest-neighbor interaction, suggesting 
the stability of dx2-y2 superconductivity against longer 
ranged interactions. 

Our findings have possible relevance for understanding 
some aspects of the effect of material specific changes in 
the electronic structure on the transition temperatures 
of unconventional superconductors. Although we have 
shown that unconventional superconductivity from repul- 
sive interactions is a robust phenomenon, which survives 
in the presence of substantial farther range repulsions, 
there is none-the-less a strong tendency for such inter- 
actions to produce a significant reduction of Tc. Con- 
versely, if screening by a proximate polarizablc medium 
reduces V and V, this could lead to a marked enhance- 
ment of Tc. 

As a point of comparison, recall that in a conven- 
tional electron-phonon superconductor, Tc depends on 
the electron-electron interactions only through ^* . How- 
ever, because of retardation, /x* ^ 1/ ln[£^i?/a;o] is 
largely independent of the bare electron-electron inter- 
action. Moreover, the effective attraction induced by the 
phonons is typically highly local, and so is also unlikely 
to be very sensitive to small changes in the environment. 
Thus, a sensitivity to the screening effects of a polarizable 
environment may be one of the uniquely characteristic 
features of pairing due to electron-electron interaction. 
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